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The existence of self-organized criticality in the Barkhausen effect and its analogy with sandpile 
models is investigated. It is demonstrated that a model recently introduced to describe the dynamics 
of a domain wall [Cizeau et al, Phys. Rev. Lett. 79, 4669 (1997)] belongs to the universality class 
of undirected Abelian sandpile models. In this way it is shown that the Barhausen effect can be 
taken as an experimental observation of self-organized critical phenomena. 
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One of the fundamental tasks of the theory of criti- 
cal phenomena is to determine the different universal- 
ity classes. If different systems or phenomena cannot be 
grouped into a reduced group of universality classes then 
the central idea of critical phenomena, the existence of 
universal behavior, will be of no relevance. This task has 
been carry out succesfully in ordinary critical phenom- 
ena. However, the precise identification of the universal- 
ity classes in a non-equilibrium critical phenomena like 
self-organized criticality (SOC) is still unresolved. 

SOC was introduced to explain the critical behavior of 
a vast class of driven dissipative system which evolve into 
a critical state 0. In its early state it was believed that 
such a critical state is insensitive to changes in control 
parameters and no fine-tuning is needed. More recent 
interpretations of this phenomena have shown that criti- 
cality in SOC systems is obtained after some control pa- 
rameters, for instance the driving and dissipation rates, 
are fine-tuned to zero Q . 

The Barkhausen effect has been taken as an experimen- 
tal observation of SOC behavior. Based on phenomeno- 
logical analogies between the Barkhausen effect and sand- 
pile models, such as the existence of power law distribu- 
tions of avalanche size and duration, some authors have 
claim that the Barkhausen effect exhibits SOC behav- 
ior |5|4|]. However, this conclusion has been criticized by 
other researches which pointed out that the observation 
of power law distributions is not necessarely an evidence 
of SOC behavior There are alternative approaches, 

like the random field Ising model , where the power 
law distributions are a consequence of the scaling prop- 
erties of disorder. The determination of the universality 
classes in this case thus become more difficult because it 
is still not clear if the Barkhausen effect exhibits or not 
SOC behavior. 

Some light in this controversy has been given by 
Cizeau, Zapperi, Durin and Stanley (CZDE) JhJ. They 
introduced an equation of motion for a single domain 
wall where dipolar interactions, demagnetization effects 
and quenched disorder are considered, containing some 
previous works as limiting cases. CZDE observed 



that their model has certain analogies with sandpile mod- 
els. More precisely the critical state is obtained when the 
magnetic field rate (the driving rate in sandpile models) 
and the demagnetization factor (the dissipation rate in 
sandpile models) are fine-tunned to zero. However their 
analysis, in relation to this analogy, was limited to this 
phenomenological observation based on numerical simu- 
lations and mean field analysis. The great importance of 
this conjecture is that if one could map the Barkhausen 
effect into certain class of sandpile models then one could 
be sure about the existence of SOC in this phenomena. 
This is precisely the scope of this work. 

We investigate the CZDE model in the case of strong 
magnetization, where dipolar interactions are relevant 
|]l3| . First we show that when the magnetization field in- 
creases at constant rate the domain wall is never pinned 
but moves in average at constant velocity. Then we ob- 
tain exact expressions for the average interface velocity 
and susceptibility. The scaling exponents are obtained 
using perturvative analysis, some of them but not all re- 
sults identical to those obtained by CZDE for the case of 
constant magnetic field. Further RG calculations reveals 
that perturbation theory is exact up to 0{r~ l ), where r 
is proportional to the square of the saturation magnetiza- 
tion. From the comparison of the scaling exponents with 
those observed in sandpile models we conclude that the 
CZDE model with magnetic field increasing at constant 
rate belongs to the universality class of directed Abelian 
sandpile models (DASM) @. 

To start our analysis let us introduce the CZDE model. 
The domain wall is modeled by a d-dimensional interface, 
dividing two regions of opposite magnetization, moving 
in a d + 1 environment described by its position h(x,t). 
Considering the contribution of magnetostatic, ferromag- 
netic and magneto-chrystalline interactions one obtains 
the following equation of motion jl0| 
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where A is a viscosity coefficient, T is the surface tension 
of the wall, H is the magnetic field intensity, L is the 
linear size of the system, and M s is the saturation mag- 
netization per unit volume. Long-range demagnetization 
effects are described by the fourth term in the right hand 
side, where Af is the demagnetization factor. Dipolar in- 
teractions are characterized by the fifth term, where the 
kernel K (x) is anisotropic and has Fourier transform 
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where 9 is the angle between k and the magnetization. 
r\ (x, h) is a Gaussian uncorrelated noise due to lattice 
defects or other factors, with zero mean and noise corre- 
lator 



(r)(x, h)r,{x', h')) = S d (x ~ x')A(h - ti), 



(3) 



where A(h) is a monotonically decreasing even function. 

If dipolar interactions and demagnetization effects are 
neglected and H is constant then eq. (0) is reduced to the 
Edwards- Wilkinson equation with quenched noise. This 
limiting case has been extensively studied in the litera- 
ture Jl5| , ^6| . A depinning transition takes place at certain 
critical field H c , determined by the disorder. For H < H c 
the interface is pinned after certain finite time while for 
H > H c it moves with finite average velocity. The up- 
per critical dimension is d c — 4. This features remains if 
dipolar interactions are considered but the upper critical 
dimension is reduced to d c = 2 . 

When the demagnetization field is included and the 
magnetic field increases at rate c then the interface is 
never pinned by impurities, but always moves with a fi- 
nite average velocity v. A perturbative solution of eq. 
(|I|) can thus be found expanding h(x, t) around the flat 
co-moving interface vt. Taking h{0, t) = vt + y(x, t), with 
{y(x, t)) = 0, we obtain the following equation for y(x, t) 
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A— y(x, t) = TV 2 y(x, t) + 2/z M s (c - 2MM s v)t 
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The average velocity is obtained using the constraint 
(y(x,t)) — 0. For this purpose is better to work with 
the equation for the Fourier transform of h(x,t), h(k,u>). 
The effective external field (c—2AfM s v)t gives a singular 
term of the order of lu~ 2 . This singular term predomi- 
nates in the low frequency limit resulting, after imposing 
(y(fe,cj)) = 0, 

(5) 



This result is valid to all orders of perturbation expan- 
sion and, therefore, exact. In the MF theory by CZDE, 
which is equivalent to the ABBM model, it is obtained 
that v ~ c/J\fM s in agreement with eq. (||). We have 
thus shown that this result is exact and, therefore, valid 
beyond the MF theory. 

Another exact result can be obtained if one com- 
putes the low-frequency and long- wavelength susceptibil- 
ity. Adding a source term (f(x, t) to the right hand side 
of eq. (Q) and going to the Fourier space one obtains the 
generalized response function 
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where 

[Goik,^)}- 1 = Tk 2 + iXuj + 4/i A/'M s 2 <5(fc, L) + K{k). 

(7) 

is the bare correlator and T,(k,ui) is the "self-energy". 
5(k,L) is the Fourier transform of L~ d . In the thermo- 
dynamic limit (L — > oo) 5(k,L) w 1 for k — > and zero 
otherwise. Since £(0,0) = and Go^O)" 1 = 4fi AfM 2 
it results that the low-frequency and long-wavelength sus- 
ceptibility (or simply the susceptibility) is given by 



X = 6(0,0) 



(8) 



This result is also exact to all orders of perturbation ex- 
pansion. 

In eq. (|8|) one cannot determine precisely which is the 
control parameter of the model. The susceptibility may 
diverges if both M or M s goes to zero. To answer this 
question we look for a self-similar solution of eq. (|J). 
Performing the scale transformation x — > bx, L — * bL, 
t — * b z t and y — > b^y where z and £ are the dynamic 
and roughness exponent, respectively, and taking into 
account eq. (|[), eq. (|J) becomes 

\b l - z ?-y{x, t) = Tb-^yix, t) - Xb^v 
at 

+b 1 ~ C r][bx, vb z t + b c y(x, t)] 
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Notice that we cannot obtain scale invariance in the elas- 
tic and dipolar intercation terms simultaneously. Thus 
we have to determine which of this terms gives the ma- 
jor contribution. A magnitude which characterizes the 
ration between elastic and dipolar interactions is 
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where A, the momentum cutoff, is of the order the inverse 
of the domain wall thickness. The case r <C 1 where 
dipolar interactions can be neglected has been already 
analyzed in |l7[ , here we focuses our attention in the op- 
posite case r> 1 where dipolar interactions are relevant. 
In this case, imposing scale invariance we obtain that TV 
and v should goes to zero and 
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The scaling exponents obtained in this way are identical 
to those derived by CZDE Jl(| for the case of constant 
magnetic field. Moreover the upper critical dimension 
d c = 2 is also the same. 

The saturation magnetization is not a control param- 
eter because the dipolar interaction term, which is pro- 
portional to M s 2 , does not renormalize. Thus, TV and v 
are the only control parameters. Since scale invariance 
is obtained when TV = + from eq. (g) one can define 
the correlation length £ ~ M~ v with v = 1. Moreover, 
from eq. (||) we obtain that the susceptibility scale as 
X ~ TV~ 7 with 7 = 1. At the critical state £ ~ L and 
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On the other hand, the increase of v may change the 
character of the noise correlator. For large v the noise 
correlator can be approximated by an annealed noise. In 
this case the magnetic field predominates over disorder, 
corresponding with a supercritical regime. From eq. (|^) 
one can define the characteristic velocity v c ~ e e with 



(13) 



which divides the phase diagram (TV, v) into two regions, 
the supercritical state v 3> v c and the subcritical one 
v <C v c . Criticality is obtained when TV — > and v — » 0. 

To go further we have performed a RG analysis of the 
problem. We integrate out the degrees of freedom in a 
momentum shell near the cutoff A and rescale k — » b~ l k, 
lo — ► b~ z tu, and w — » b^ +d+z w, where b — e l with I — > 0. 
The flow equations for the parameters T, A, TV, M s and 
v are obtained through a direct application of the RG 
transformations to eq. (Q) in the Fourier space. The 
renormalization of the moments of the noise correlator 
(Qn = J A(q)q n ) is obtained considering vertex func- 
tions [ ^6|l^ ]. As in perturbation theory here we only 
consider the case r » 1. In this case we found correc- 
tions of the order of r -1 to the exponents computed by 
perturbation theory p8| . Thus, the exponents z and Q in 
eq. ( O ) are exact up to 0(r _1 ). 

Let us now determine the avalanche distribution ex- 
ponents. The Barkhausen signal V(t) is the voltage pro- 
duced from a pickup coil around a ferromagnet subjected 



to a slowly varying applied field. In the low-frequency 
limit the time scale for domain wall motion is much 
smaller than the time between jumps and, therefore, one 
may guarantee that each induced voltage jump corre- 
sponds with a single avalanche in the domain wall mo- 
tion. A resolution voltage level Vr is defined, such that 
one can not resolve details below Vr. An elementary 
Barkhausen jump can thus be defined as the portion of 
the V(t) signal delimited by two subsequent intersections 
of the signal with the Vr line. With this definition, the 
duration T is simply the time interval between these two 
subsequent intersections and the size s is the area delim- 
ited by V(t) and Vr between the same points. 

In the subcritical regime the dynamics takes place in 
the form of avalanches, characterized by the avalanche 
size P(s) — s~ T f(s/s c ) and duration P{T) — 
T~ a g(T /T c ) distributions, where s c and T c are the 
avalanche size and duration cutoffs. In the subcritical 
state s c ~ (T x l a and £ ~ t~ v while at criticality s c ~ L D 
and £ ~ L, where D is the avalanche dimension, leading 
to the scaling relation 
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Other scaling relations are obtained taking into account 
that x = (s) and / dsP(s) = J dTP{T), which lead to 

7 = ^— ^, ( T -l)D = (a-l)z, (15) 
a 

respectively. 

Using these scaling relations and our result 7 = 1/ = 
z = 1 we obtain the following scaling relations for the 
avalanche exponents 
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D, 



(16) 



On the other hand, for d < d c , the avalanche dimension 
and the roughness exponent are related via D = d + £, 
while above the upper critical dimension one obtains 
D = d c = 2. Hence, we can compute r and a using 
eq. (K&j and this value of D. For instance 



a = I , for d 
= 2, for d ■ 



1 = 2; 
1 > 3. 



(17) 



The scaling law in eq. (^2|), the scaling relations in 
eq. (116 ) and the values of the scaling exponent a in 
eq. (]l7) are also obtained for DASM jl4|]. However, 
the upper critical dimension of DASM is 3, and not 2 
as obtained here. This apparent contradiction can be 
understood if one takes into account that in DASM time 
evolution does not introduce an additional dimension be- 
cause it can be represented by the evolution in the pref- 
erential dimension. On the contrary in the motion of the 
domain wall the time evolution take place perpendicular 
to the (i-dimensional substrate and, therefore, introduces 



3 



an additional dimension. Taking this fact into account 
the upper critical dimension will be 2 + 1, as in DASM. 
Thus, the CZDE model in d + 1 dimensions is equivalent 
to a d-dimensional DASM. 

The fact that the CZDE is mapped into a directed 
sandpile model, and not into an undirected one, is due 
to dipolar interactions which introduces an anisotropy 
in the system. If dipolar interactions becomes negligi- 
ble (r <C 1) then we expect that the CZDE model will be 
mapped into an undirected sandpile model. Experiments 
in magnetostrictive materials, where dipolar interactions 
can be neglected, have been performed (l^ jT^,Q . Earlier 
measurements by Urbach et al gives r = 1.33 ± 0.10. 
More recently, Durin and Zapperi (DZ) |2^] reported 
the more accurate exponents r = 1.28 ± 0.02 and a = 
1.5 ± 0.1. On the other hand, numerical simulations of 
the CZDE model in two dimensions and without dipo- 
lar interactions have also been performed Jl^,|2(|. The 
more accurate numerical estimates, reported by Durin 
and Zapperi @, are r = 1.26±0.04 and a = 1.40±0.05. 
These exponents are in the range reported for undirected 
Abelian sandpile models in two dimensions. For instance, 
numerical simulations of the Bak-Tang-Wiesenfeld (the 
prototype of undirected Abelian sandpile model) gives 
t = 1.293 and a = 1.480 @. However, to validate this 
guess we cannot base our analysis only in the avalanche 
exponents, we must provide further elements. In Jl7t we 
have shown that the CZDE without dipolar interactions 
can be mapped into a class of undirected sandpile models 
with annealed noise in the toppling rule. 

In conclusion we have shown that the CZDE model for 
the dynamics of a domain wall belongs to the universality 
class of directed Abelian sandpile models. Dipolar inter- 
actions introduce an anisotropy in the systems leading 
to a preferential direction along the magnetization. We 
have provided strong arguments which states the equiv- 
alence between the dynamics of a domain wall and sand- 
pile models. This conclusion can be extended to systems 
with many domain walls, since for short time scales the 
interaction between domain walls can be neglected and 
the single domain wall picture is correct. Hence, the 
Barkhausen effect actually exhibits self-organized criti- 
cal behavior. 
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